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Phase Synchronization with Type-II Intermittency in Chaotic Oscillators 
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We study the phase synchronization (PS) with type-II intermittency showing ±2tt irregular phase 
jumping behavior before the PS transition occurs in a system of two coupled hyperchaotic Rossler 
oscillators. The behavior is understood as a stochastic hopping of an overdamped particle in a 
potential which has 2-7r-periodic minima. We characterize it as type-II intermittency with external 
noise through the return map analysis. In et < e < e c (where et is the bifurcation point of type-II 
intermittency and e c is the PS transition point in coupling strength parameter space), the average 
length of the time interval between two successive jumps follows (l) ~ exp(| e t — e | 2 ), which agrees 
well with the scaling law obtained from the Fokker-Planck equation. 
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Synchronization is one of the basic phenomena ubiqui- 
tously found in physical, chemical, biological and physi- 
ological systems. In the classical sense, synchronization 
of 'periodic self-sustained oscillators is usually defined as 
locking of the phases, nQ\ — mQ 2 = const, due to weak 
interaction while the amplitudes can be quite different. 
This phenomenon has been quite well studied and has 
witnessed a lot of practical applications in various engi- 
neering fields 

Recently, the notion of synchronization has been ex- 
tended to coupled chaotic oscillators (i.e., individual os- 
cillators are chaotic without coupling). One of the re- 
markable developments is the observation of phase syn- 
chronization (PS) phenomenon in a system of two mu- 
tually coupled nonidentical self-sustained chaotic oscilla- 
tors . The phenomenon is analogous to synchroniza- 
tion of periodic oscillators where only the phase locking is 
a matter of importance. Above a critical strength of cou- 
pling, suitably defined phases of two chaotic oscillators 
lock each other and synchronize, while their amplitudes 
remain chaotic and uncorrelated with each other. Also 
it was found that the phase difference 4> between two os- 
cillators increases with an intermittent sequence of 2ir 
jumps before the PS transition occurs. It means that 
phase slips occur from time to time and the phase dif- 
ference changes by 2tt during a rather small interval of 
time. 

The intermittent behavior and its scaling properties 
near the PS transition in a coupled Rossler system was 
studied by several authors They provided an ex- 

planation for the phase jumps by reducing the original 
system into a simplified model describing an overdamped 
particle sliding in a "noisy wash-board potential" . Also 
by studying the scaling rules of the jumping behavior, 
they found that the phenomenon is related with type-I 
intermittency M in the presence of noise Hpl- So far, 



this has been known to be the only route to the PS tran- 
sition in two coupled self-sustained chaotic oscillators, 
while nothing forbids other types of intermittency to ex- 
ist. 

In this paper, we report another route to PS transi- 
tion exhibiting ± 2n jumps, which is characterized by 
the type-II intermittency with external noise. We con- 
sider the following two coupled hyperchaotic Rossler os- 
cillators (HRO's), 

xi.2 = -^1,22/1,2 - zi,2 + e(x 2 ,i - Xi j2 ), 

Vl,2 =^1,221,2 + 0.252/1,2 + 101,2, 

Zl.2 = 3.0 + 211,231,2, (1) 

wi.2 = —0.5^1,2 + 0.05wi,2 + e(w 2 ,i - 101,2), 

where two variables x and w are coupled and the sub- 
scripts 1 and 2 refer to each of the oscillators. Here 
^1.2 = 1.0+AO/2 is the overall frequency of each chaotic 
oscillator, and e is the coupling strength. 

To observe PS we must define a suitable phase re- 
lated to this system. Since the phase portrait of a hy- 
perchaotic Rossler oscillator in the x-y plane explicitly 
shows a rotational trajectory around a center, (xo,yo) 
as shown in Fig. 1, the phase can be defined in relation 
with this rotation, i.e., 2:1,2 and j/1,2 are transformed into 
polar variables A\ 2 and Si, 2 around the center of the 



rotation where A 1:2 — \J (2:1,2 — xq) 2 + (3/1,2 



yo) 



and 
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yi.2-yo 



v XI, 2 — XQ 

can be rewritten as follows 



With these variables, Eq. (1) 



Ai,2 = 0.25^4i,2 sin 2 #1,2 - 21,2 cos 6>i, 2 + wi,2 ' 
+ e(A 2 .i cos 6*2,1 cos 6*1,2 - ^1,2 cos 2 6*2,1), 
01,2 = ^1,2 + 0.25 sin 6*1,2 cos 6*1,2 
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COS 02,1 — COS^i^ ) Sin 01,2, 



(2) 
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z l,2 = 3.0 + (j4i )2 008 01,2 + x )zi t 2, 

"i,2 = -0.5zi, 2 + 0.05wi, 2 + e(iV2,i - u>i,a)- 

where 101,2 is a shorthand for wi,2 + £0^1.2 and xq = 
— 16.0 and y — 0.0 are used. By solving Eq. (1) or (2) 
numerically, the phase difference (cf) = 81 — 2 ) of the two 
oscillators is obtained as shown in Fig. 2. Note that the 
onset of PS occurs near e c = 0.18 ||. Before the onset 
of PS, the smaller the coupling is the more frequently 2-7T 
jumps occur like the case of the Rossler system. But in 
our system, the phase difference jumps by 2tt upwards 
or downwards irregularly as shown in Fig. 2, while in 
the Rossler case it jumps only upwards monotonously. 
In this respect, our system should be distinguished from 
the Rossler system qualitatively 0,[|J|. 

Here, we are interested in this new type of phase jump- 
ing behavior and would like to uncover the basic struc- 
ture behind it. For this purpose, the following equation 
for the phase difference cf> is obtained from Eq. (2): 



= AO + (0.25 + e) (sin 0i cos0i -sin0 2 cos0 2 ) 

+ — — sin 0i — — sin 02 + -7— cos 0i — — cos 02 
A\ A2 A\ A2 

— e ( — ^ sin 0i cos 02 j- cos 0i sin 02 J . (3) 

\A\ A 2 J 

The r.h.s. of Eq. (3) can be rewritten in terms of <p by 
grouping relevant terms, and then we obtain the following 
phase equation similar to that describing phase locking 
of periodic oscillators in the presence of noise |2|,|],[§ : 



= Afi + a sin <j> + (5 sin ~ + £, 



(4) 



where 



a = (O.25 + e )cos(0 1+ 02)-|^ + ^]^ 
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Note that there are two time scales, i.e., the fast one, 
27r/f2i,2, related with the frequency of each individual 
oscillator in Eq. (2), and the slow one, 2n/A£l, originat- 
ing from the frequency mismatch in Eq. (3) which is the 
characteristic time of <f> dynamics. Since a and /3 and £ 
are fast fluctuating pieces compared with the time scale 
of 4> variable, qualitative features of 4> dynamics can be 
revealed after averaging a and (3 over the slow time scale, 
while £ is left intact as an external noise in Eq. (4): 



(•5) 



where 



F(4>,e) = An + (a) sin, 



</?>sin|. 



(6) 



In Eq. (5) and (6), by <f> we mean a new phase vari- 
able which simulates original (f> dynamics. Here (a), ((3) 
represent averaged values for a long enough time inter- 
val which is order of 27r/Ar2. This equation describes 
an overdamped particle moving in a potential under the 
influence of external noise £, where the potential V(<fi, e) 

is defined by F(4>) = §0. This averaged po- 

tential is obtained by integrating the force F((j>,e) with 

respect to 0, i.e., V((f>, e) = — J F(<fi', e)d(j>' up to an 
arbitrary integration constant: 



V{4>, e) = -AtlcP + (a) cos $ + 2</3) cos |. 



(7) 



Figure 3 (a), (b), and (c) show time series of <f> at 
e = 0.16, the corresponding force F, and the potential 
V respectively. As shown in Fig. 3 (c), the existence 
of 27r-periodic minima is the main feature of the poten- 
tial related with the phase jumping dynamics. Note that 
the overall slope of this potential is negligible due to the 
small constant force term Af2, so that the "wash-board 
potential" picture Q is not applicable to studying the bi- 
furcation related with the intermittency in this system. 
With this potential structure, however, we can qualita- 
tively explain the phase jumping dynamics as a hopping 
of an overdamped particle between 27r-periodic minima 
in the potential through the stochastic process driven by 
the external noisy force £. In Fig. 3 (c), when a particle 
is in a potential well, it can jump randomly to the near- 
est neighboring potential wells since the barrier heights 
of both sides of the well are similar. When the well is not 
deep, the particle resides there for a short time interval, 
but it resides there longer when the well is deep. This 
behavior can be easily observed in the time series of <j). 

To determine the type of intermittency explicitly, we 
obtain the return map of <fi from the original Eq. (2) 
by finding <j> n at every 27r rotation of phase 0i. Fig- 
ure 4 (a) shows a schematic relation among force, poten- 
tial, and the corresponding return map. In Fig. 4 (b), 
the return map shows many cells separated by 2tt corre- 
sponding to 27r-periodic minima in the potential struc- 
ture. The trajectory in the return map moves from one 
cell to the neighboring cells irregurally in accordance with 
±27r jumps. To see the structure more clearly, we ob- 
tain the return map of (j> from the averaged Eq. (5) in 
one cell at every 2tt rotation of phase 0i when £ = 0. 
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The solid curves in Fig. 4 (b), (c), and (d) are the re- 
sults and they quite well follow the one obtained directly 
from the original system (2). This means that the aver- 
aged Eq. (5) properly describes the main features of ±2n 
phase jumping dynamics. From the shape of the curves, 
it is clear that the return map is actually that of type-II 
intermittency ]7|Jli|]. In Fig. 4 (c) the center is a stable 
fixed point which becomes unstable when the slope at 
this point equals 1.0 at the bifurcation point e t = 0.038 
as shown in Fig. 4 (d). The trajectory can escape from 
the center to the nearest neighboring stable fixed points 
located each at ± 2n distance away through the stochas- 
tic process. Therefore the scaling of the average laminar 
length of 4> (here the laminar length is the time elapsed 
between two successive ±2ir jumps) is expected to follow 
that of type-II intermittency in the presence of external 
noise. 

In order to verify the above conjecture, we compare 
the average laminar length scaling relation obtained nu- 
merically from Eq. (2) with the analytical result ob- 
tained from the Fokker-Planck equation (FPE). The local 
Poincare map of type-II intermittency with external noise 
is described by the following difference equation [jl0| |l3| : 

x n+1 = (l + e)x n + axl + V2D£ n , (8) 

where a is a positive arbitrary constant, and D is the 
dispersion of Gaussian noise £ n . It is well known that 
the above map (8) can be approximated to the backward 
FPE in the long laminar region fl4| : 

9G(x,t) , , dG(x,t) n 9 2 G(x,t) 

dt ~ 1 ' dx + 8 2 x ( ' 

where G(x,t) is the probability density of particle at 
(x, t). The scaling relation for the average escaping time 
can be derived analytically from this FPE. According to 
the analytical estimation made in our recent work p5| , 
the average laminar length scales as follows when e t < e: 

(I) ~ (Ooexp(|e t -e| 2 ). (10) 

Now, we perform numerical simulation to obtain the 
scaling law of the average laminar length from Eq. (2) 
and the results are shown in Fig. 5. First we get the 
value of y-intercept, ln(Z)o, in Fig. 5 (a) to determine the 
scaling exponent by the linear regression. The slope of 
the regression line in Fig. 5 (b) is very close to 2.0 within 
± 5 % error, which agrees with the scaling exponent in 
Eq. (10). Therefore our conjecture based on the return 
map analysis is verified and strengthened by this remark- 
able agreement with the result obtained from the FPE. 
Also this agreement in turn justifies our averaging out 
the fast fluctuating factors a and (3 in Eq. (5), which 
marginally contribute as multiplicative perturbations in 
phase jumping dynamics. 

In conclusion, we have firstly found the type-II inter- 
mittency route to the PS transition in a system of two 



coupled self-sustained HRO's. In this system, the phase 
difference between two oscillators exhibits apparently ir- 
regular ± 27r jumps near the onset of the PS transition. 
Furthermore, we have identified this novel behavior as 
type-II intermittency with external noise through the re- 
turn map analysis and the FPE approach. 
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FIG. 1. The phase portrait of one of two coupled HRO's 
projected on the x-y plane. Initial values for HRO's 1 and 
2 are xi(0) = -20.0, yi(0) = zi(0) = 0.0, Wi(0) = 15.0 
and x 2 (0) = -20.1, y 2 (0) = z a (0) = 0.0, w 2 (0) = 15.1 
respectively. Eq. (1) is numerically solved by using a 
fourth-order Runge-Kutta method with the parameter mis- 
match Afl = 0.001, and the same values are used throughout 
this Letter. 
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FIG. 3. (a) Time series of phase jumping dynamics, (b) 
the time-averaged force F ((/)), and (c) the potential V((f>) at 
e = 0.16. 




e = 0.160 . 








, 










:'. (c) 









e = 0.160 














(b) . 


-2% 




e t = 0.038 ,;' 






.'• -(d) 



FIG. 4. (a) schematic relation among force, potential, and 
return map; (b) return map (black dots) of original system 
(1) and return map (gray dots) from the averaged Eq. (4); (c) 
and (d) magnified return maps around one potential well with 
e = 0.16 and et = 0.038 respectively. Note that the tangent 
line at the center stable fixed point becomes tangential to the 
diagonal line in (d). 



FIG. 2. Time series of phase difference <f> in two coupled 
HRO's for various values of (a) e = 0.16, (b) e = 0.17, and 
(c) e = 0.18. ± 2-7T phase jumps are shown for the cases of (a) 
and (b), while for (c) the PS transition seems to be achieved 
for the observed time scale. 
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FIG. 5. Scaling relation: average laminar length \n(l) vs 
\e t - e | 2 in (a) and ln(ln(Z) - 2.6) vs. In | e t - e | in (b). Note 
that the slope in (b) is approximately 2.0 within ± 5 % error 
in agreement with the analytic estimation from the FPE. 
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